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ABSTRACT 


Many current army combat simulation models require 
digitized terrain inputs for their execution, resulting 
in a costly and time-consuming operation. A methodology 
to generate parametric terrain as a function of the | 
number of hill masses and the height, location, spread, 
and slope of the hills is presented in this paper. Modi- 
fied bivariate normal (MBVN) distribution function and a 
linear function cutting the MBVN function are employed to 
achieve the terrain representation. Analysis of various 
terrain parameter values was conducted and the significant 
trends resulting from the study are described. Finally, 
recommendations for future enhancement of the model are 


given. 
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I. INTRODUCTION 


Many current army combat simulation models require 
digitized terrain inputs for their execution. This process 
is very costly and requires a significant amount of computer 
Storage. In addition, computation of line-of-sight and 
mobility require a great deal of computer running time. 

The current methodology utilized to represent terrain 
is also questionable from a statistical point of view. 

A Vector Research Incorporated (VRI) study on terrain line 
Of sight, concluded that "... present and past army study 
results, based on the analysis of combat results on a very 
limited sample of terrains, may have been determined by 
the terrain selection process and not by the actual weapon 
system or force design differences." More specifically 

it stated the following: 

A. "There is extreme sensitivity in combat model 
results as the scenarios (terrain and movement assumption) 
are varied, even when variation is within a class of 
scenarios chosen for their a priori equivalence." 

B. "This sensitivity can be slightly reduced, but 
remains extreme (with probabilities of win estimable only 
within plus or minus 25%) even when battle results are 


1 
used to redesign scenarios." 


tearrell, Robert L. and Richard J. Freedman, Investigation 
of the Variation of Combat Model Predictions with Terrain Line 
Sifsaght, Vector Research Incorporated, 1975. 
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These results imply that sufficient replications of each 
type of terrain should be run in order to reach a satisfactory 
statistical level of significance. For VRI's analysis at 
least fifty replications of each type of terrain were used. 

In view of the problems stated above, research was under- 
taken by Major Christopher Needels* to develop a new method- 
Ology for the representation of terrain. In particular, a 
modified bivariate normal distribution function was utilized 
to generate hill masses in a wide variety of configurations. 
Line of sight and movement routines were also developed to 
evaluate various intervisibility parameters over many terrain 
configurations for very low cost. 

This thesis represents an extension of the models 
developed by Needels to generalize the terrain generation 
routines. In particular, methodology was developed to 
generalize the hill mass generation from symmetric hill 
shapes to non-symmetric configurations. The methodology 
developed for this generalization is described in Chapter III. 

In order to evaluate the effects of the number of hills, 
hill height, and slope on various intervisibility parameters, 
an experiment was designed to measure these effects. A 
description of the design methodology and results of the 
analysis are given in Chapter IV. Finally, the conclusions 
drawn from the current analysis and recommendations for 
future research are presented in Chapter V. 


-Needels, Christopher J., Parameterization of Terrain in 


Army Combat Analysis, Master's Thesis, Naval Postgraduate 
School, Monterey, California, 1976. 





If. REPRESENTATION OF TERRAIN 


A. HIGH RESOLUTION COMBAT MODEL REPRESENTATION 

The classic approach to representing terrain in land 
combat models is to utilize digitized terrain evaluation 
data from Waterway Experiment Station (WES) for the 
particular area of interest. 

The Dynamic Tactical Simulator (DYNTACS) model serves 
aS a good example of the state of the art in terrain 
representation and use. DYNTACS is a two-sided dynamic 
model of battalion level combat representing the detailed 
interaction of elements in a combined arms environment. 

DYNTACS utilizes 100 meter grid squares for which 
macro-terrain elevations are determined from WES tapes. 
The simulation divides each grid square diagonally, thus 
providing a series of adjoining triangular terrain. The 
entire battlefield, therefore, is represented as a surface 
of equal size triangles which vary in slope depending on 
the elevation at their corners. 

Needels (pp. 9-10) provides a concise description of 
the DYNTACS methodology for line of sight calculation as 
follows. 

From the macro-terrain data, line of sight 
between any two opposing elements is computed. 
This is accomplished by first computing the 
angle between the horizontal and a straight 
line drawn between an observer and target 
(O-T Line). The program then conducts a 


search of the terrain along the path of the 
O-T line to see if any macro-terrain is 





higher than the O-T line itself. This is 
accomplished by comparing the angle of the 
O-T line with the angle above horizontal 

of the Observer-Terrain line. If the latter 
angle is larger, there is no intervisibility. 
Over the duration of a battle with numerous 
elements this calculation may be made 
thousands of times. Consequently, not only 
is the time to prepare the tapes high, but 
also the time to compute lines of sight once 
the terrain is input to the model. 
Considerable effort by developers and users 
of this model has been expended to stream- 
line this subroutine. 


B. SYMMETRIC REPRESENTATION 

The basic models developed by Needels* (hereafter referred 
to as SIMTER) which provide the basic foundation for the 
methodology developed in this thesis are described in this 
section. The basic motivation for SIMTER was to develop a 
random terrain representation which could be replicated 
quickly but at the same time be representative of a particular 
“type" of terrain. The approach adopted for SIMTER generated 
a variable number of hill masses. The shape and height of 
the hills could be varied, but each hill was symmetric, since 
it was generated directly from the modified bivariate normal 
(MBVN) density function. 

SIMTER has two parts. The first part is the main 
program which creates terrain and, at Eve discretion of the 
user, plots both a three-dimensional drawing and a contour 
map. Used strictly for terrain generation, it can be used 
as a preprocessor for other combat models by producing grid 


points and elevations similar to those provided on computer 


3Needels, Parameterization of Terrain in Army Combat 


Analysis. 





tapes by WES. The second part of SIMTER is a movement/line 
of sight routine which moves a target along specified routes 


across the generated terrain. 


C. A NEW REPRESENTATION METHODOLOGY 

The methodology developed in this thesis provides for 
the generation of hill masses utilizing MBVN density 
function to represent masses that are not symmetric. This 
development provides a significant increase in flexibility 
to represent hill masses more realistically. The mathematical 
development of this methodology is presented in the next 
chapter, after which actual results from executing the model 


are given. 


ak 








IIil. THE EXTENDED MODEL 


A. UNSYMMETRIC TERRAIN 


The unsymmetric slope of terrain can be generated by 


a cutting method using the MBVN as follows: 








L X-U xX“ oynH ve 
Z = MH-EXP } - x a)? - 2 BR) t (7 
m0) x x Yy iv, 
where MH : maximum height of hill 


Ze hex@ght of hail 
u nu? spread of hill 


a8 ellipse factor of hill 


The unsymmetric slope of the hill can be represented by the 
vertical value of the MBVN subtracted from the vertical value 


of a linear function as shown in Figure 1 by the value HI. 





FIGURE 1. UNSYMMETRIC SLOPE GENERATION 
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Two basic questions must be addressed: 
1. Where do we start the cutting of the MBVN function? 
2. How do we compute the initial cutting position? 
If Mis the direction of the cutting line (Line A), then 
the cutting plane is through the location of the center of 


mae Hill (Wy Hy) as shown in Figure 2. 







er ! 


(XIN, YIN) yy 


FIGURE 2. CUTTING LINE DETERMINATION 


femene Initial point of cutting is (XIN,YIN) and HI is 

height at the initial point, this point is on the cutting 
plane and also on the MBVN distribution function. Thus, if 

we consider the graph in two dimensions as shown in Figure 2, 
Line B is tangent to the circle at the point that the vertical 
value on the MBVN function 1S HI and the slope of Line B 
tce—iyte Lf the initial point o£ the cutting is (XIN,YIN) 
this point is on Line B, Lame A and the circle at a height 

HI on the MBVN function. We can compute the initial point 


(XIN,YIN) as follows. 


1S 





im ‘o@ = 0 Case 
XIN-y YIN-y 
HI = MH*EXP(- = (—_——#) 7 + 14 
Oo O 
_ Y 
XIN—-y Veen = 
-2 In(HI/MH) = ( X24 2 
OF oy 
From (2), the derivative is computed as follows: 
AIN—yp 
Semanal 
ee _ Ox ee 
aXIN YIN-w M 
ee 
av 
From (2) and (3) 
XIN = y+ o, 2 
1L+M (5/5) 
-—2 1 I /MH 
YIN = u, #0 M 2h 
(o,./9,,) +M 
ne ee O., then 
N= uw. ¢ o., -2 In(HI/MH) 
- 1+M 
= MH 
Age SEG eae eee, 
Y Y 1+M 
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(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


i) 





If the hill is cut from left to right and oe Dl then 





-2 1ln(HI/MH) 
1 +M* 







XIN = yp -o (8) 






-2 ln(HI/MH) 


(9) 
1 + MC 


YIN = i =-o M 
y 


imewe hill is cut from right to left and Tareas OR ss then 







= 2 ln HME) 





XIN = u +o (10) 
a as 1+ M- 

ey = Lu + 0M 72 in(HI/MH) Gab) 
Y Y 1+ M- 





Thus, if the direction of cutting and HI are selected, then 
the starting point can be computed. Care must be taken in 
the selection of the value of HI. If HI is too small, then 
the hill will have an unrealistic shape. The value of HI 


affects the shape of the hill as shown in Figure 3. 








. wn 
eel 
| 7, ‘os 
co } aie 
7 = aie iy pes ot 
/ ae —_ ~ ae a 
| . E y re 
} Le / Es oN 
aos d oa \ 
| 7 mn Me H =. \ 
Se | ee oie 
—e ea Pi ae Y 
| ee ae a" X 
Aa ” ~ 
ae i ze | i | een 7 
x e ~ = 4 4 ee: - te 5 ' na 


FEGURE 3. EFFECT OF HI ON THE HILL SHAPE 
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Ze Smo JeeCace 


XIN-y XIN-p. YIN-wu 
HI = MH * EXP{ - —+-[(——¥)?-2 (—_—¥) (__¥) 
2(1=-p~) x ox Cy 
YIN-u 5 
+ (———*) G2) 
O 
y 
XIN-y XIN“. - YIN-u 
ieee a yin(kt/Ma) = (———%)? — 2 (—_—_*) (___) 
O O GO 
x x y 
YIN-wu 5 
+ ( 5 ) (13) 
v4 
From (13) the derivative is computed as follows: 
1 XIN-w, YIN- 
a aa ian = )] 
x x y 
a4IN _ ee 3 
dXIN — ~ M ee 
AIN-w YENa i 
- = [-o———*) + (1 
x ¥ 
From (14) 
XIN- o. + Moo oN 41 
(as ee (—-—*) (a5) 
O 
x Mo,, + po, 4 
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Feom (13), (15) and (16), 









An = 
-2(1 ae ln (HI/MH) = (1 = 20 a +) X) 2 
E 2 = 
C x 
Let DB = l1- 20 + i 
C 2 
C 
Then 
mIN =< uito | r2gbee%) in (Ht/oy 
uw, = Oy . 
Also, let 
a 2 
ie - C + 20C + l 
Then, 
~2(1- 0%) 1n(HI/MH) 


E 
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(18) 


be!) 


(20) 


(21) 





FIGURE 4. CUTTING DIRECTION 


If the cutting direction goes through the longer 


axis as shown in Figure 4, then 


2 Z 
Oy =F 0 
2pc ow 


The cutting direction, a, is determined by 


‘ a 2 7 a7 
Q = +s arccot 
y 


2 =2 00 16 (23) 
Xx 


If the cutting direction goes through the shorter 


axis as shown in Figure 4, then 


8B = 90° -a (24) 
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Next, the height of intersecting hill masses must 


be determined as illustrated in Figure 5. 





FIGURE 5. INTERSECTING HILL MASSES 


HB defines height of hill "B" after cutting. 
HA defines height of hill "A" after cutting. 

Height of hill "A" is greater than that of hill "B" 
at the point XA but height of hill "B" is greater than that 
of hill "A" after cutting. If there is more than one hill 
located close to each other, those hills must be adjusted 
by comparing the height after cutting each one. The 
computation of the height of the terrain at any point is 


described in the following section. 


B. TERRAIN HEIGHT 

The height of the hill which is generated by the MBVN 
distribution corresponds to the height of a hill which is 
generated by the cutting method. The height of a hill 


at any point is computed as follows: 


IO, 








The height of the MBVN density function is given by 


Cah 





x- X-W you yY-u 
2 = MHEXP {- —=—— [(——%)? - 29(—%y(—_¥) +) 2} 
= ej 0 ej 
2(l-p”") ae x y sé 
(25) 
The height of the linear function is 
ZCUT = (X-XIN)tan 6+ (Y - YIN) tan 9 (26) 


@D 
I 


where cutting angle of the x-axis 


cutting angle of the y-axis. 


oe 
II 


The height of the hill after cutting is 


0 if 22 ana - 2ZCUTREcrO 
0 if Z°< HE and. ZCUT > 6 
Z-HI shag Zi re eeecin Clem ZO Uden 


Z- ZCUT - HI ii Tr An ZOU. => =O 
and Z > ZCUT+HI 


0 1£ 2 SHE and ZCUT > 0 
and Z < ZCUT+HI 


If more than one hill exists at a point on the terrain, the 


greatest height after cutting is used. 


Ce eemevEMENT LINE OF SIGHT 
The procedure of computing the line-of-sight between 


any two points on the terrain is described in this section. 


20 





If the position of the observer is (OX,OY) and the position 


of the target is (X,Y), then the slope of the Observer- 
Target (O-T) Line is given by 


(28) 


If the intersection point of a hill and the O-T Line is 


(XI,YI), then HI and the O-T Line are perpendicular as 
shown in Figure 6. 





oY} —~ 
f 
: 
\ 
FIGURE 6. INTERSECTION OF O-T LINE 
WITH HILL MASS 
Tas , . 
ML? ©OX -ML -OY+uU_t+ML-u 
x 
XI = 5 (29) 
1 + ML 
Vue =  ME(XT = OX) + OY (30) 
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If there is more than one hill on the O-T Line, the 
dominant hill (after cutting) is used to compute the percent 
of the target covered relative to the observer. Target 
movement along the surface slope of the hill is computed 


as follows: 





ny ot! 
aH = MH*EXP {~ —"—>— [(——*)* ~ 29(- >) (Se 
2(l-p -) Xx x 4 
ou 
+ (=X) 71} 
Y 
- (X-XIN) tan 6 - (Y-YIN) tan 6 (31) 
92H ~MH (oy p(y~Hy) 
= 1-92 a Toy 








X~ XxX“ yru y-u 
*exp{——+ | (<8)? -29 (— 8) (=) + (2) 71) ~ tan 8 

















0 0 
2(1- x y 
Gee) x ae 
92H ~MH ss p(y-Hy) 
a «(| ( 
oy l-o¢ oe Oo y 
, = xe yal Vallee 
*exp (2 (8) * -20( 9) (4) + (1) tan 
Z(l=0°) ‘ x y 


iz 





If DEL is the direction of target movement, then 


= 02H 92H ; 
ae COs x + a Sate ve (34) 


If the speed of the target is constant "V" on the surface, 
then 


velocity = us 


eee (35) 
J 1 + (DEL) * | 


If the speed of the target is dependent on the slope of 


the surface, then 
velocity = V + W(DEL) 


{36) 


where W is the speed factor of slope. 
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IV. TEST AND EVALUATION OF MODEL RESULTS 


A. DESIGN AND CUTTING ANGLE METHODOLOGY 

The methodology for parametrically describing terrain 
as a function of the number of hills, the center and spread 
of each hill, and the slope was developed in the previous 
chapters. In this chapter the effect of the number of hills, 
height of hills, and the slope on various intervisibility 
measures is investigated. Three random terrains were 
generated (with each terrain having six, twelve, and eighteen 
hills) resulting in nine basic terrain configurations. 

For each configuration three average hill heights (100, 
200, and 300 meters) were utilized. 

Finally, for each of the 27 resulting configurations, 
four slope cutting methods were used as follows: 

1. Down slope-cut from observer side to target side. 

2. Side slope~-cut from right to left as viewed from 
the target. 

3. Up slope — cut from target side to observer side. 

4. Both slope-symmetric Slope as viewed from any 
position. 

A total of 108 unique terrain configurations were 
generated for the investigation. For each Gf the cenrigd— 
urations, the mean and standard deviation describing the 
center and spread of each hill is specified. In addition, 
the cutting angle for determining slope must also be 


specified. 
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If the cutting angle is constant, then the hill can be 


described as follows: 

HILL = F(Peak, Position, Standard Deviation, Slope) 
In the constant cutting angle case, only the slope is 
variable, with peak, standard deviation, and position 
constant. In other words, if the slope varies in Figure 7, 


the standard deviation of the resulting distribution also 


varies. 


A LINE 


FIGURE 7. CUTTING ANGLE AND SLOPE 


The objective of the methodology is to maintain a constant 


standard deviation of the distribuition function . resulting 
from each cutting angle, so that only the slope of the hill 
varies. Otherwise the determination of how slope affects 

the intervisibility parameters would not be possible, since 


it would be confound with the standard deviation (spread) 


Guetene hills. 
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The problem, therefore, was to determine a method to 
use various cutting angles and maintain a constant standard 
deviation. After extensive investigations into various 
approaches to this problem, no method was found which would 
mathematically guarantee a constant standard deviation. 

The method shown in Figure 8 was adapted as the best 


approximation method available. 





FIGURE 8. CUTTING ANGLE METHODOLOGY 


In Figure 8, HT is the maximum height of the hill after 
cutting with Line A. Note that the cutting angle is now 
specified in reference to the center line of the original 
bivariate normal hill mass distribution. In order to 
generate a symmetric hill mass with approximately the same 
standard deviation (Line C) the cut is made maintaining 
a constant maximum height, HT, as in Figure 8. 

Several test runs made using this methodology indicated 
that the resulting standard deviations were very nearly 
equal for cutting angles. However, very little change in 


slope results for cutting angles greater than twenty degrees. 
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B. OBSERVER-TARGET POSITIONS 

For each terrain an observer position and initial target 
position were selected. Two routes from the target to the 
observer position were selected in accordance with accepted 
tactics of advance routes. Selected terrains indicating 
the observer-target position and selected routes are shown 
in Figures 12 through 17, In order to assure that the 
"number of hills" factor was properly analyzed, the same 
routes of approach on each generated terrain were used, 


independent of the number of hills. 


C. TARGET MOVEMENT 
Two cases of target movement were considered in the 
analysis. 
fee constant Speed 
The target moves along the route at a constant speed 
of six meters per second, independent of hill slopes. The 
constant speed runs serve as the base case for measurement 
of the percent of time intervisible along the route. 
2. Variable Speed 
The speed of the target is a function of the terrain 
slope being traversed. The speed of the vehicle is computed 


by the relationship as follows: 
Velocity, si— 6.0 5 4000s lope) 


This slope is the directional derivative value at the 


instantaneous target position on the hill. 
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D. TARGET INTERVISIBILITY 

Five dependent variables describing target intervisibility 
are computed for. each advance route as follows: 

1. Percent of time intervisible along the route, 

2. The number of intervisibility segments along the route, 

3. The average intervisibility segment length, 

4. The average intervisibility segment time, and 

5. The total length of intervisibility over the route. 
Intervisibility parameters are computed by utilizing the 
line of sight algorithm at fixed time intervals, TL, during 
the movement run. For this analysis, TL was set equal to 
one second. A sample output for one compute movement trace 


is given in Figure 9. 


E. EXPERIMENTAL DESIGN 

The objective of the analysis is to determine the effects 
Semene mumber Of hills (A), the height of"hilis (B) and the 
slope of the hills (R) on target intervisibility. To 
accomplish this objective, an experiment was designed as 


described by the following design matrix. 


SLOPE 
+ of Hills Height Down Side Up Both 


12 200 


18 200 
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EEGURE 9. 


OUTPUT OF DATA 


Definition of Variables 


Number of hills 
Spread of hills 
Maximum height of hill 


Different amount of the spread and height 
among hills 


Grid increment for the terrain matrix 
Size of the grid system in meters 
Center, pesiticon of hai 


Maximum height of hill in modified bivariate 
normal distribution 


Spread of hill (standard deviation) 
Cutting angle 

Cutting direction 

Initial position of cutting plane 


Hieght of hill at the center position on MBVN 
after cutting 


Coordinates of approaching route of target 


Movement time 


MAP COORDINATES Target position at given time 


ELEVATION 


SPEED 


VISIBLE . 


HLIN 


Target elevation on the terrain 

Velocity on the terrain at the given time 

1 means visible and 0 means invisible 

The different height of O-T line at the 
maximum height of surface of the terrain on 


the O-T line, and the lower height between 
observer and target 





DZ 


KL 


DLIN 


Maximum height of surface of the terrain on the 
O-T line 


The hill which is maximum height between the 
observer and target 


Difference maximum height of surface between 


observer and target, and lower height between 
observer and target. 
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The appropriate linear model is as follows: 


: oe = +A. +B. +R + .. + 
ijk eee ee ee ABR, TT 


+ Fiske 


“ash = intervisibility dependent variable measure 


u = grand mean of the population 
A = number of hills, where 
6 hades 
1 = 2-12 hills 
3-18 hills 
B = average height of hills, where 
1-100 meters 
J = 2-200 meters 
3-300 meters 
R = direction of slope of hills, where 


l-down slope 
2-side slope 
3-up slope 
4-both (symmetric) 
Within each cell of the matrix, there are six replications 
consisting of two routes on each of the three basic terrain 
configurations. 

For the first treatment levels of slope defined above, 
a cutting angle of fifteen degrees was used, the symmetric 
case was generated with a cutting angle of zero degrees 
using the methodology previously described. Four analyses 
were conducted as representative of the large number of 


possible sitatuations. The measurements selected were 


1. Percent of time intervisible for constant speed 
and for variable speed, and 


2. Number of intervisibility segments for constant speed 
and for variable speed. 
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The results for each selected analysis are described in 


the following section. 


F. ANALYSIS OF RESULTS 


1. Percentage of ENtervisibziaey Time 


a. Constant Speed 


os Transgeneration Data) 
SUMMARY TABLE - ANOVA 
Source DF SS MS F Ratio 
Ort hills 2 0.0437 0.0218 2366417 
A 
height of hill 2 G20 ay, 0.0058 6.3209 
B 
steppe of hill 3 0387 O27 029 bs) SOE 
R 
AB interaction 4 0.0006 0.0081 O2bs27 
AR interaction 6 0.0026 0.0004 0.4759 
BR interaction 6 0.0016 0.00027 0.2943 
ABR interaction 12 0.0012 0.0001 Oe Fisk 
Beror 179 O-1655 0.0009 


Each of the three main effects above are significant 
at the 0.01 level, but no interaction terms are significant. 


The following table gives the mean and standard deviation 


for each treatment level of slope. 
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Slope Mean Std. Dev. 


down 
Side 
up 

both 


0.6130 
0.5214 
Une 26 


2573 


0.1539 
Oey cal 
0.1629 
0.1595 


Note that the up slope has the greatest percentage of inter- 
visibility time; down slope, both, and side slope following 
iietaat Order. 

The following table gives the mean and standard 
deviation for the treatment levels of number of hills and 


height of hills for the constant speed case. 


# of Hills Mean Standard Deviation 
6 0.6976 0.1620 
12 0.5462 O76 20 
18 Qv5272 0.1490 
Height of Hill Mean Standard Deviation 
100 0.6475 Oe 72 
200 0.5844 0.1684 
300 O7.55.5.6 Oj.6:2.5 


35 





b. Speed which Depends on Slope 


3/8 


(x Transgeneration Data) 


SUMMARY TABLE - ANOVA 


Source DF SS MS 


F Ratio 
fOr nilis 2 Or 5152 0.2576 34.286 
A 
Heigimt of hill] 2 0.2041 0.1020 13.562 
B 
Slope of hill 3 0.3830 OL 277 16.991 
R 
AB interaction 4 0.0033 0.0068 0.1090 
AR interaction 6 0.0381 0.0064 0.8450 
BR interaction 6 0.0165 0.0028 0.3670 
ABR interaction 12 0.0091 0.0008 0.1000 
Error 179 14 seks 0.0075 


Each of the three main effects above are significant at the 


0.01 level, but no interaction terms are significant. 


c. Constant and Dependent Speed. 
The ANOVA for comparing constant and dependent 


speed is given in Figure 10. The treatment level data is 


as follows: 


Mean Standard Deviation 
Constant Speed 0.569 Oy ulay pe 
Dependent Speed ORs) Shs 6: Calo Zl 
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2. Number of Intervisibility Segments 


a. Constant Speed 


SUMMARY TABLE - 


(xl /2 Transgeneration Data) 


Source DF 

%# Of hill 2 
A 

Hemet Of hill 2 
B 

Slope of hill | 
R 

AB interaction 4 

AR interaction 6 


ABR interaction eZ 


Peror 179 


Sys) 


8.3415 


0.2402 


Phos PEAT Ee: 


0.0659 
0.2686 
Cel oe9 


9.0475 


ANOVA 


MS 


4.1708 


One 20u 


02.9095 


0.0165 
0.0448 
0.0116 


0.5056 


F-Ratio 


82.281 
Bo, ES, 
IL ees 


8) 5 SaaS: 
0.8838 
0.2288 


From the ANOVA, the number of hills and slope are significant 


at the 0.01 level, but the height of hills is not significant. 





b. Dependent Speed 


- SUMMARY TABLE - ANOVA 


iE 
(X 2 Transgeneration Data) 


Source DF SS MS F-Ratio 

# of hill 2 Ba 2224 Ns il Mi big 19.337 
A 

Hemant Of hill 2 0.2561 0.1281 ZeAZ72 
B 

Stepe of hill 3 2.6384 0.8795 16.972 
R 

AB interaction 4 0.06446 0.01611 Ors VE 

AR interaction 6 0.2905 0.0484 0.934 

ABR interaction 2 0.1289 0.0107 0.2065 

Error 179 9.2756 0.0518 0.2065 


The results of the dependent speed case are analogous to 
those of the constant speed case. That is, the number of 
hills and slope of hills are significant, but the height of 
Rells is not significant. 

The mean and standard deviation of the treatment levels 


are given below. 


oi) 








SLOPE MEAN STD DEV 
CONST DEP CONST DEP 
DOWN W207 5 ee On 1.04444 1.0444 
SIDE ARG LAG 1.1339 1.1339 
UP Bp SUL cl choles 1.1947 Ie Alou 
BOTH 4.4259 4.4259 1.3818 1.3818 
# OF HILL MEAN STD DEV 
CONST DEP CONST DEP 
6 4.930 3.038 eS 7 1 0.7166 
ee 4.139 4.139 0139 1.0388 
18 3.042 4.931 Om 2 0r7 i alg feu: 
HEIGHT OF MEAN STD DEV 
ae CONST DEP CONST DEP 
100 3.875 Be 1.393 1.393 
200 4.069 4.069 1.191 1b aL). 
300 are 5 4.169 1.191 1.195 
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V. CONCLUSION AND RECOMMENDATION 


Terrain is one of the key factors in the analysis of 
land combat. Thus, the representation of terrain in a land 
combat model is a critical factor in the model design. High 
resolution combat models typically employ actual terrain 
data to describe the terrain environment. This approach, 
although realistic for the specific terrain being represented, 
1s very costly and time consuming to implement. 

The goal of the research by Needels and this thesis is 
to provide a parametric representation of terrain usable 
either as input to a high resolution combat model or as a 
"stand-alone" model to investigate intervisibility segment 
length parameters. In particular, the current model 
describes the terrain as a function of the number cf hill 
masses, the height of the hills, the position and spread 
of the hills, and the slope. 

A MBVN distribution function and a linear function 
cutting the MBVN function are employed to achieve the 
terrain representation. The terrains given in Figures 12 
through 17 illustrate the flexibility of the model. In 
fact, with appropriate selection of the input parameters 
(given in Figure 9) any desired terrain may be generated. 

The results of the analysis described in Chapter IV 
indicate that slope is highly significant in the resulting 


values of intervisibility parameters. Some of the trends 
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noted from the results given in Chapter IV are described 
below. 

1. An up slope (cut from target side to observer side) 
resulted in the largest value of percent intervisibility 
time. 

2. Increases in the number and height of hills resulted 
in decreases in the percent intervisibility time. 

3. Speed dependent on terrain slope resulted in a 
slight reduction in percent intervisibility time from the 
constant speed case. 

The analysis of the hypothetical data described in 
Chapter IV serves to illustrate the potential utilization 
of the model. The trends indicated from the analysis are 
appropriate for the selected terrains and conditions, but 
should not be construed as generalized results. 

An alternate method (tilted method) for generating 
unsymmetric terrain outlined below is suggested as an area 
for future enhancement of the model. Given a MBVN density 
function centered at Y in Figure 11, define 6 as the tilted 
angle, assuming 9 = 0. The height of the density function 


16 22 at Y2 and is given by 





X-U yru 
22 = MH*EXP{- > r(——%) 2 + (—— 4%) 2773 (38) 
oO ) 
x Y 
and 
ee tan eo mer), (39) 
Y2-Y 
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FIGURE lL TILTED METHOD 


From 38 and 39, 


oe You 
MH*EXP (- 5 [(——%)* + (—¥P = (y-y,) tan(Z- a) (40) 
x 


ve 





For a specified value of x and y, the value of Y> can be 
computed using the golden section method. The tilted height, 


ZH, can then be computed as follows: 


LES. 


Sin 0 (41) 


ZH 





Additional investigation of this method is required before 
implementation in the model. 

In conclusion, a model has been developed and exercised 
which provides a potentially valuable tool for efficient 
analysis of various terrain effects in the land combat 


environment. The intervisibility parameters output from 


the model can serve as critical inputs to land combat models 
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for evaluating such things as the contribution of mobility, 


agibility, and firepower to battlefield survivability. 
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